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We study the thermodynamic properties of a dirty ferromagnetic S|F|S Josephson junction with s-wave su- 
perconducting leads in the low-temperature regime. We employ a full numerical solution with a set of realistic 
parameters and boundary conditions, considering both a uniform and non-uniform exchange field in the form 
of a Bloch domain wall ferromagnetic layer. The influence of spin-active interfaces is incorporated via a mi- 
croscopic approach. We mainly focus on how the electronic specific heat and density of states (DOS) of such 
a system is affected by the proximity effect, which may be tuned via the superconducting phase difference. Our 
main result is that it is possible to strongly modify the electronic specific heat of the system by changing the 
phase difference between the two superconducting leads from up to nearly n at low temperatures. An en- 
hancement of the specific heat will occur for small values h c± A of the exchange field, while for large values of 
h the specific heat is suppressed by increasing the phase difference between the superconducting leads. These 
results are all explained in terms of the proximity- altered DOS in the ferromagnetic region, and we discuss 
possible methods for experimental detection of the predicted effect. 

PACS numbers: 85.25.Dq,74.25.Bt,74.45.+c,74.78.Na 



I. INTRODUCTION 

In recent years, due to the important role of hybrid struc- 
tures with superconducting and magnetic layers in vital 
circuit-elements like transistors and high-resolution devices 
like detectors, such structures has attracted much attention 
from the research community. In this way, several interesting 
phenomena of such systems have been predicted both in the 
dirty and clean limits of transport, and subsequently been ob- 
served in experiments: non-monotonic dependence of the crit- 
ical tempe rature T r in S|F hybrid structures on the thic kness 
of F layer ' 1 ' 2 ' 3 ' 4 ' 5 ' the 0-tt transitions in S|F|S junctions,™! 
and the appearance of odd-frequency pairing correlations^^ 
just to mention a few. The main cause of the appearance 
of these interesting phenomena is the proximity effect be- 
tween the superconductor and the ferromagnet, where Cooper 
pairs leak from the superconducting side to the ferromagnetic 
layer. S|F|S junctions are routinely fabricated by experimen- 
talists these days, and currently such hybrid structures are 
intensely investigated due to their potential both in terms of 
functionality 11 and novel fundamental physics that may be 
exploredp2E3 

In the context of applications, studies of the thermody- 
namic properties of superconductors have mostly focused on 
electron cooling propertied, although recently the influence 
of the proximity effect on the entropy production in a non- 
magnetic Josephson junction was investigated 15 . A key to un- 
derstanding the thermodynamic properties of a system is the 
behavior of the DOS near Fermi level, and any control param- 
eter that can adjust the DOS in an efficient and well-defined 
manner would offer significant advantages with respect to tai- 
loring desired thermodynamic properties. 

Very recently, it has been studied numerically 16 and demon- 
strated experimentally 17 how the density of states (DOS) may 
be altered controllably in such structures by creating a non- 
magnetic S|N|S Josephson junction and generating a super- 



current. For a ferromagnetic Josephson junction, however, it 
remains to be clarified precisely how the DOS is influenced 
by the phase difference in the full proximity effect regime. A 
new feature which is expected to come into play for ferromag- 
netic Josephson junctions is the presence of odd- frequency 
superconducting correlations, which can induce a qualitative 
shift in the DOS from a low-energy minigap- structure^ to an 
enhancement 19 . 

In this paper, we show how it is possible to obtain a huge 
enhancement of the specific heat of a ferromagnetic Joseph- 
son junction at low temperatures, simply by tuning the su- 
perconducting phase difference by means of either a current 
or an external magnetic flux in a SQUID-like geometry. We 
demonstrate explicitly how the predicted effect occurs for a 
set of realistic experimental parameters, and how it persists 
even in the presence of an inhomogeneous magnetization tex- 
ture such as a Bloch domain wall in the ferromagnet. We find 
that the enhancement of the specific heat is strongest for ex- 
change fields h comparable in magnitude to the superconduct- 
ing gap A, i.e. h ~ A, whereas for higher exchange fields the 
effect eventually vanishes. Our findings can be verified exper- 
imentally by using calorimetry techniques or high-resolution 
thermometry,^ and could have potential applications in de- 
vices utilizing an active tuning of the thermodynamic prop- 
erties of nanoscale conductors. We underline that while it is 
well-known that the DOS in a Josephson junction is sensitive 
to the phase difference, our main result pertains to the manner 
in which the DOS varies and the resulting consequences for 
the electronic specific heat of the junction. 



II. THEORY 

To investigate the physical properties of the S|F|S Joseph- 
son junction, one alternative is to solve the quasiclassical 
Eilenberger equation^! to obtain the Green's functions. In 
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FIG. 1: (Color online) The S|F|S junction with i) uniform ferro- 
magnet in the z direction and ii) with a Bloch domain wall in the 
ferromagnet. The magnetization texture for the Neel wall would be 
obtained by replacing the x -component of the magnetization with an 
^/-component in case ii). The red arrows show the magnetic moments 
in the F layer. 



the diffusive limit, the Eilenberger equation reduces to a sim- 
pler set of equations known as the Usadel equations 26 . For 
numerical purposes, it is convenient to use a parameterization 
method for the Green's functions. One approach for the pa- 
rameterization is the (9-parametrization 23 as follows: 

A _ / M cog + (M • a)s p+ \ 

g ~\ p~ -Mocao-(M-ays)' 

= c[i(M z 02 - Myas) ± M x ao] ± M ois, (1) 

where Oj are the identity (j = 0) and Pauli (j = 1,2,3) 
matrices, and 

a = (01,02,03). (2) 

Also, s = sinh(0) and c = cosh(0). The Green's function is 
then completely determined by the complex functions 6, Mo, 
and M with the additional constraint Mq — M 2 = 1 in order 
to satisfy g 2 = 1. The other approach for a parameterization 
of the Green's functions is the Ricatti-parameterization 24 25 . 
We found that for our purposes in this paper, i.e. a full nu- 
merical investigation of the density of states and consequently 
the thermodynamic properties of a diffusive S|F|S junctions, 
this parametrization is much more numerically stable than the 
^-parameterization. The Ricatti-parameterization can read as 
follow™ 

/A/;(l-77) 2AT 2 \ 

^"V K{-i+rd)' 

Here, g 2 = 1 since 

M= (l^jf)- 1 M= (I + 77)" 1 . (4) 

We use for 2 x 2 matrices and . ? . for 4 x 4 matrices. In 
order to calculate the Green's function g, we need to solve 
the Usadel equation^ with appropriate boundary conditions 
at x = —dp j 2 and x = dp/2. We introduce the supercon- 
ducting coherence length as £s = -s/DsJAq. Following the 
notation of Ref. 27 , the Usadel equation reads 



and we employ the following realistic boundary conditions for 
all our computations in this paperP^ 

2(d F gdg = [g B cs (</>), 9] + i(Gs/G T )[diag(73, 73), g] (6) 

at x = —dp/2. Here, d = and we defined ( = Rb/Rf as 
the ratio between the resistance of the barrier region and the 
resistance in the ferromagnetic film. The barrier conductance 
is given by Gt, whereas the parameter Gs describes the spin- 
dependent interfacial phase-shifts (spin-DIPS) taking place at 
the F side of the interface where the magnetization is assumed 
to lie in the yz-plane, being parallel to the z-axis at the inter- 
faces. The boundary condition at x = dp/ 2 is obtained by 

letting G s -> (~G S ) and # B csO) -> [-#bcs(-</>)] in Eq. 
([6]), where 

7 BCS (<A)= + c)e'*/ 2 , 

W*)=W*) e_1 *- C7) 

Above, is allowed to be different from Gs in general. 
For instance, if the exchange field has opposite direction at 
the two interfaces due to the presence of a domain wall, one 
finds Gs = —Gs- The total superconducting phase differ- 
ence is (j), and we have defined s = sinh(#),c = cosh(#) 
with $ = atanh(A /£^) using A as the superconducting gap. 
Note that we use the bulk solution in the superconducting re- 
gion, which is a good approximation when assuming that the 
superconducting region is much less disordered than the ferro- 
magnet and when the interface transparency is small, as con- 
sidered here. Effectively, the inverse proximity effect is thus 
ignored. We use units such that H = ks = 1. 

The values of Gs and Gt may be calculated explicitly from 
a microscopic model, which allows one to characterize the 
transmission {t J n a } and reflection amplitudes {r J n a } on the 
j G {5, F} side. Under the assumption of tunnel contacts 
and a weak ferromagnet, one obtains with a Dirac-like barrier 
modePSEl 

Gt = GqY^Tu, G s = 2G Q £ [p F n - ^) (8) 

n n 71 

upon defining T n = ^ a |^ a | 2 , pi = Im{r£ T (r£ a )*} and 

= Im{£^ t(^n,|)*} anc * also for simplicity, we assume that 
the interface is characterized by N identical scattering chan- 
nels and consequently omit the subscript 'n' in the summa- 
tions. 

To find the specific heat of the system, we need to calculate 
the local density of states normalized against its normal-state 
value 

N(x, E, T, 4>) = Tr(Re[Af(l - jj)})/2. (9) 

We assume that the S electrodes are not influenced by the 
proximity effect, the total electronic specific heat (C tot ) of the 
SFS junction can be determined by 



Dd(gdg) + i[Ep 3 + dmg[h a,(h- af},g] = 0, (5) 



C tot (T,d>) = C F (T,<p) + C s (T). 



(10) 
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(a) d p /^ s =0.3 and h/A^O.5 



(b) d p /^ s =0.3 and h/A^IO 




(c) d F /^ s =1.0 and h/^0.5 





E/A, 

(d) d p /^ s =1.0 and h/^slO 




FIG. 2: (Color online) The normalized local density of states vs. energy and phase difference between two s-wave superconducting leads for 
a uniform magnetization texture. We have presented several values of the thickness of ferromagnetic layer <1f/^s and exchange field h/Ao, 
and the DOS is evaluated in the middle of the junction, i.e. x = 0. 



Here, Cs(T) is the specific heat of superconducting plates 
while 



C F (T^)=TdS F (T^)/dT, 



(11) 



is the specific heat of the ferromagnetic part of the junction. 
The entropy of the ferromagnet layer in the proximity system 
can be obtained from 



"(4/i) J 



L/2 poo 

dx / dEN(x,E,T,(j))x 

L/2 JO 



(12) 



{f(E) ln[f(E)] + [1 -f(E)] ln[l -/(£)]}, 

and f(E) = {1 + exp[£/T]} _1 is the Fermi-Dirac quasipar- 
ticle distribution function at temperature T. 

Since we employ a numerical solution, we have access to 
study the full proximity effect regime and also, in principle, an 
arbitrary spatial modulation h = h(x) of the exchange field. 
This is desirable in order to clarify effects associated with non- 
uniform ferromagnets, such as the presence of Bloch domain 
walls. 



In this paper, we will consider two different types of mag- 
netization textures: homogeneous magnetization and a Bloch 
domain wall structure as shown in the Fig. [TJj and ii), re- 
spectively. It is seen from Fig. [TJpart ii) that for a domain 
wall structure, the magnetic moment has two components un- 
like the homogeneous type. The Bloch model is given by 
h = h(cos6y + sinOz) and its structure is shown in Fig. 1 
part ii), where we defined = — 2 arctan(x/<i^) 25 , with dw 
as the width of domain wall. Moreover, the center of the F 
layer is located at the origin x = 0. Below, we shall consider 
a domain wall of width dw I dp = 0.5, thus ensuring that the 
magnetization is fully aligned with the z-axis at the interfaces. 



III. RESULTS AND DISCUSSION 

In this section, we present our main results of the paper, 
namely the manner in which the DOS of a S|F|S diffusive 
Josephson junction is altered in the presence of phase differ- 
ence between two s-wave superconducting leads and conse- 
quently the behavior of the electronic specific heat of such 
junctions. We consider a ferromagnet with two types of mag- 
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(a) d p /^ s =0.3 and (b) d p /^ s =0.3 and h/A Q =1 




FIG. 3: (Color online)The normalized local density of states vs. energy and phase difference between two s-wave superconducting leads and 
inhomogeneous magnetization texture for several values of thickness of ferromagnetic layer cLf/^s and exchange field h/Ao. The DOS is 
evaluated in the middle of the junction, i.e. x = 0. 



Homogeneous exchange field Domain wall 




FIG. 4: (Color online) The normalized electronic specific heat of the diffusive S|F|S josephson junction vs. phase difference between two 
s-wave supperconducting leads with a ferromagnetic layer featuring a homogeneous (left panel) and inhomogeneous Bloch domain wall 
magnetization texture (right panel). The arrows indicate an increasing exchange field h. 
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netization texture: homogenous and inhomogeneous (Bloch 
domain walls), including also the role of spin-active inter- 
faces at the two boundaries. In the quasiclassical framework 
employed here, we have to consider an exchange field much 
weaker than the Fermi energy in order to remain within the 
regime of validity. For a weak, diffusive ferromagnetic al- 
loy such as PdNi, the exchange field h/A can be varied from 
a few meV to tens of meV by changing the relative content 
of Pd and Ni. Even weaker exchange fields h of order meV 
are found in for instance Y4C03, Y9C07, and TiBei.gCuo.J^. 
Therefore, we shall here consider exchange fields ranging 
from 0.5 meV up to 5 meV. The scenario of a thin junction 
will be contrasted with that of a thick junction 
d>F/ f £s = 1-0. For a superconducting lead like Nb with a 
coherence length of £s — 18 nm > the ratio of <ip/£s=0.3 pro- 
vides a ferromagnetic layer thickness equal to 6 nm which is 
experimentally accessible 32 . The temperature will be fixed at 
T=0.05T C and consequently our results are valid for low tem- 
perature regime. The spin-dependent interfacial phase-shift 
(spin-DIPS) term — Gs / Gt is obtained via the micro- 
scopic theory introduced in the previous section, and depends 
e.g. on the magnitude of the exchange field and the interface 
transparency. Since we calculate Gs microscopically for a 
simplified model with a Dirac tunneling barrier, Gs is not 
treated as a phenomenological parameter here. We choose 
[ip = 1 eV and us = 10 eV for the Fermi level in the ferro- 
magnet and superconductor, respectively, and consider a rela- 
tively low transparency barrier of Z = 3. The electron mass 
rriF and in both of the F and S regions is taken to be the 
bare one (^0.5 MeV). The ratio of the electronic resistances 
of the barrier region and the ferromagnet layer is assumed to 
be (=Rb/Rf=^ throughout our computations. We also in- 
sert a small imaginary part S=10~ 3 x Ao into quasiparticle en- 
ergies i.e E — > E + iS for access to more stability in our 
computations. The small imaginary part can be interpreted 
as accounting for inelastic scattering. As we discuss below, 
we find that the specific heat of the S|F|S diffusive junction 
can be strongly enhanced by changing the phase difference 
between the two singlet superconducting leads from up to 
values near tt for both a homogeneous exchange field and in 
the domain wall case. Due to limitations of our numerical 
code, we were not able to investigate phase differences (j) very 
close to 7r. The huge enhancement of the specific heat can 
be seen even for values of the exchange field several times 
the superconducting gap in the domain wall case. Upon in- 
creasing the magnitude of the exchange field further to values 
h ^> Ao, this effect vanishes. The enhancement is most re- 
silient towards an increase in h in the case where a domain 
wall is present. Both the enhancement of the specific heat and 
its persistence in the domain wall case can be understood by 
investigating the DOS in the ferromagnetic region. We now 
proceed to a presentation of our main results. 



A. Density of states (DOS) of a S | F | S junction at low 
temperatures 

In this section, we discuss the behavior of the DOS in a fer- 
romagnetic region by changing the phase difference between 
two superconducting leads connected to it. We fix the temper- 
ature at T=0.05T C and also use from microscopically values 
for spin-DIPS term in the two boundaries. The results are 
shown in Fig. [2] for the homogeneous exchange field case, 
while the domain wall scenario is demonstrated in Fig. [3] 
In both figures, we provide a contour-plot of the DOS in the 
middle of the F layer as a function of quasiparticle energy 
E measured from Fermi level and the superconducting phase 
difference (j). 

Let us first consider the homogeneous case shown in Fig. 
ga) for the case ^f/^5=0.3 and ft/A=0.5. The most obvious 
feature is that a minigap- structure is induced in the low-energy 
regime close to the Fermi level, flanked by a peak structure 
below the gap and at the gap. The minigap is shown to close 
as the phase difference moves towards <p = ir, as is also the 
case for S|N|S junction^. In Fig. (2^c), the junction thick- 
ness is increased to dp /is = 1-0, and it is seen that the peak 
structures remain. The main difference from (a) is that the 
low-energy DOS is enhanced, indicating the odd-frequency 
correlations are present and comparable in magnitude to the 
even-frequency correlations. The minigap is split into two 
and is seen to shift away from zero energy. The appearance 
of the multiple peak structures as a function of energy E orig- 
inates from an effective superconducting gap felt by each spin 
species which is different in magnitude for spin- j and spin- 
| quasiparticles. This is similar to the scenario of thin-film 
superconductors subjected to an in-plane external magnetic 
field 33 . When the exchange field is increased to Ji/Aq = 10 as 
shown in (b) and (d), only the standard BCS-coherence peaks 
remain at the gap, although much weaker in magnitude than 
in the bulk case. 

We now turn to the domain wall case, shown in Fig. [3] The 
most noteworthy change from Fig. [2] is that the zero-energy 
DOS is enhanced in (b) and (c). This observation signals that 
odd-frequency correlations are stronger in the domain-wall 
case, a finding which agrees with the results in Ref. 27 . The 
physical reason for this is that the inhomogeneous magnetiza- 
tion texture generates not only the 5^ = triplet component, 
but also the long-ranged S z = ±1 triplet components, which 
also are odd in frequency due to the isotropization caused by 
the impurity scattering. 

The presentation of the DOS and its dependence on the en- 
ergy E and superconducting phase difference (j) presented in 
this section is a useful preliminary which, as we shall see, ex- 
plains the origin behind our main result of a strongly enhanced 
specific heat, which we shall now move on to. 



B. Electronic specific heat of S|F|S junction at low 
temperatures 



Let us consider the electronic specific heat of the S|F|S dif- 
fusive Josephson junction vs. phase difference of the two 
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superconducting leads. As mentioned in the Introduction, 
the phase difference is an experimentally tunable quantity by 
means of e.g. current-biasing the junction or applying an ex- 
ternal magnetic field in a SQUID-like geometry. Since the 
proximity effect is in general much weaker in the low-energy 
regime for dp /is — 1.0, we focus here on the more interest- 
ing case <ip/£s=0.3. 

The results are shown in Fig. [4] As seen, the left panel 
is related to the homogeneous exchange field scenario while 
the right panel is related to the inhomogeneous magnetiza- 
tion in the form of a Bloch domain wall. In both cases, the 
curves show a giant enhancement of the normalized specific 
heat when the exchange field is comparable in magnitude to 
the superconducting gap. For larger exchange fields, the spe- 
cific heat becomes a monotonic, nearly constant function of 
the phase difference 0. We note that the enhancement per- 
sists for larger values of h in the domain wall case (up to 
h/Ao ~ 3.0) compared to the homogeneous case. The phys- 
ical reason behind the enhancement of the specific heat stems 
from the dependence of the DOS on </>, as shown in Figs. [2]and 
[3] For instance, for a very weak exchange field h/Ao = 0.5, 
the DOS-plots presented in the previous section showed how 
the minigap closed with increasing <j). Since it is the low- 
energy DOS that mainly contributes to the electronic specific 
heat, increasing the phase difference <p will naturally lead to 
an increase in Cp. More specifically, we have verified numer- 
ically that at T/T c =0.05, only energies up to E/A ~ 0.35 
contribute to the specific heat integral in Eq. §V2) . 

By increasing the magnitude of the exchange field in the 
ferromagnetic layer, a kink appears in the specific heat. To 
identify the cause of the appearance of the kinks, one should 
investigate the related DOS of the system. Consider now 
for concreteness the DOS in the homogeneous exchange field 
case with h/A =l.l, which is seen to display a kink in the 




4>/tt 



FIG. 5: (Color online) The normalized local density of states of 
the diffusive S|F|S junction vs. phase difference and energy a for 
homogeneous structure of ferromagnetic layer with exchange field 
h/Ao — 1.1 and thickness (If Us — 0.3. 



specific heat in the left panel of Fig. [4] The kink of this curve 
appears near (fr/ir ~ 0.7, consequently leading us to plot the 
DOS of the system near this value vs. E/A and </>/7r. The 
resulting DOS is shown in Fig. [5] As seen, the cause of ap- 
pearance of a kink near (fr/ir ~ 0.7 is the zero-energy peak 
that occurs in this region of the phase difference. Such a zero- 
energy peak should be a direct result of the manifestation of 
odd-frequency correlations in the system 19 . In Fig. [5] it is 
seen that an abrupt conversion takes place at c^/tt ~ 0.7 along 
the E = line from a fully suppressed DOS to an enhanced 
value compared to the normal- state. Such an abrupt conver- 
sion was also very recently studied in RefP^, where it was 
demonstrated that the conversion was associated with a tran- 
sition from pure even-frequency to pure odd-frequency corre- 
lations. The simultaneous decrease of the DOS when mov- 
ing away from the Fermi level results in a rapid decrease of 
the specific heat, thus leading to the non-monotonic behavior 
shown in Fig. [4] 



IV. SUMMARY 



In summary, we have considered the density of states and 
electronic specific heat of the diffusive S|F|S Josephson junc- 
tion both for a homogeneous and inhomogeneous magnetiza- 
tion texture, including the role of spin-active interfaces. We 
find that the electronic specific heat of the S|F|S junction can 
be tuned to undergo a strong enhancement by increasing the 
phase difference between two superconducting leads. The ex- 
perimental requirement for observation of this effect is that 
the width dp of the ferromagnetic interlayer is considerably 
smaller than the superconducting coherence length (typically 
dp in the range 5-10 nm), and that the exchange field is com- 
parable in magnitude to the gap. The effect persists in the 
domain wall case up to exchange fields h/A ~ 3, yielding 
h in the range 4-7 meV for a weak ferromagnetic alloy. Our 
prediction may be tested by e.g. calorimetry measurements 
of the Josephson junction, and the results reported here could 
have interesting consequences for nanoscale devices relying 
on an active tuning of their thermodynamic properties. 
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